Residual Stresses in Glasses 
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The history dependence of the glasses formed from flow-melted steady states by a sudden cessation 
of the shear rate j is studied in colloidal suspensions, by molecular dynamics simulations, and mode- 
coupling theory. In an ideal glass, stresses relax only partially, leaving behind a flnite persistent 
residual stress. For intermediate times, relaxation curves scale as a function of 'jt, even though no 
flow is present. The macroscopic stress evolution is connected to a length scale of residual liquefaction 
displayed by microscopic mean-squared displacements. The theory describes this history dependence 
of glasses sharing the same thermodynamic state variables, but differing static properties. 
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Materials are often produced by solidification from the 
melt, involving nonequilibrium quenches. This imprints 
a history-dependent microstructure that strongly affects 
macroscopic material properties. One example is resid- 
ual stresses [H 12 : if particle configurations cannot fully 
relax to equilibrium, some of the stresses, arising in the 
presence of flow in the melt, persist in the solid. 

Small glass droplets (known as Prince Rupert's drops 
or Dutch tears since the 17th century) vividly display the 
effects of residual stresses [3]: they withstand the blow 
of a hammer onto their main body, but explode when 
the slightest damage is inflicted upon their tail (releasing 
the frozen- in stress network). Today, safety glass and 
"Gorilla glass" covers for smartphones are deliberately 
pre-stressed during production to strengthen them. A 
theoretical understanding of residual stresses and their 
microscopic origins is however still not achieved. 

We seek to understand generic mechanisms by which 
residual stresses arise. A convenient starting point is to 
investigate the stress relaxation a{t) following the ces- 
sation of shear flow of rate 7, from a well-defined non- 
equilibrium stationary state (NESS). Such "mechanical 
quenches" are ubiquitous in soft matter, where pre-shear 
is applied to "rejuvenate" the otherwise ill-defined glassy 
state [IHZ]. For these systems, the soft-glassy rheology 
model (SGR) [5] predicts asymptotic power laws that im- 
ply the relaxation of stresses to zero [5] ■ In the following, 
we will reserve the term residual stress to describe a fi- 
nite, persistent stress remaining in the (ideal) glass even 
at arbitrarily large times after the cessation of flow. 

In addition to macroscopic rheology, we investigate the 
evolution of the microscopic dynamics as characterized 
by the waiting-time dependent mean-squared displace- 
ments (MSD). The latter reveal the dynamical shrinkage 



of shear-fluidized regions after cessation, and phenomena 
akin to, yet different from the intensely studied aging 
dynamics after thermal quenches [TOl [11] . 

Experiments on a variety of colloidal suspensions, to- 
gether with molecular-dynamics (MD) simulations, pro- 
vide a coherent qualitative picture that can be ratio- 
nalized by mode-coupling theory of the glass transition 
(MCT) [T2] within the integration-through-transicnts 
(ITT) formalism |13j . The theory in particular predicts 
the existence of a residual stress in the glass, the magni- 
tude of which depends on the history of shear. We expect 
the same mechanisms to be generically valid for many - 
colloidal as well as molecular - glass formers. 

The experiments, theory, and simulation are described 
in detail in the Supplementary Information (SI) [14] and 
summarized here. We perform rheology on colloidal poly- 
styrenc-PNIPAM core-shell particles in aqueous solution 
(PP) dSmS], and on PMMA particles in different sol- 
vents (HS) [TTHT^ . These are well studied, nearly hard- 
sphere glass formers |20j . but differ in particle proper- 
ties like softness and polydispersity. The average size 
is i? ~ 90 nm (PP) and R 267 nm (HS) for macro- 
scopic rheology, and R « 770 nm (HS) for confocal mi- 
croscopy. Density is expressed as a dimensionless packing 
fraction ip; the glass transition occurs in the two systems 
at ~ 0.64 (PP) and (fc « 0.59 (HS). In the PP system, 
the effective packing fraction is sensitively tuned through 
temperature. 

MD simulations using a dissipative-particle-dynamics 
thermostat were performed on a binary glass-forming 
Yukawa mixture as outlined in Refs. [191 [21]. The unit 
of length was chosen as the small-particle radius R. The 
glass transition temperature of this system is « 0.14. 

ITT-MCT was evaluated numerically [HMj with 
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FIG. 1. Stress decay a{t) after cessation of steady shear 
(time in units of the microscopic relaxation scale tq), for var- 
ious shear rates -yro — Peo (increasing from red to blue) and 
control parameters, as labeled. Dotted lines indicate times 
used to extract residual stresses in Fig. [5] (a) MD simulation: 
T = 0.14 in the liquid, T = 0.1 in the glass, (b) Isotropic 
hard-sphere model of ITT-MCT, (^mct = 0.51, 0.515 (liq- 
uid), and 0.52 (glass), (c) HS colloidal suspension: ip — 0.542, 
0.587 (hquid), and 0.614 (glass), (d) PS-PNIPAM particles: 
T ^ 18°C {(f ^ 0.57, liquid, 7 as labeled with g = 3.4) and 
T = 15 °C (v? ^ 0.65, glass, g = 4.0). 

the Percus-Yevick model of the equilibrium hard-sphere 
structure (giving ipc « 0.516) and an isotropic approx- 
imation to spatial integrals. This model (ISHSS) has 
been used together with our MD system and confocal 
microscopy to study the evolution from equilibrium to 
the NESS [ini |25] . Mean-squared displacements are cal- 
culated using a schematic model based on Refs. |27] . 

Figure [l] shows the transient decay of the shear stress 
a{t) measured in rheology and computer simulation, and 
calculated within ITT-MCT. For each system, curves for 
various shear rates 7, and thermodynamic control vari- 
ables above and below the glass transition are shown. 
Stresses are reported in entropic units, kT/R^ « 0.21 Pa 
(HS), 5.33 Pa and 5.82 Pa (PP at T = 15 °C and 18 °C). 
Values measured in the two colloidal suspensions differ 
by a factor 25 in their absolute value, consistent with 
their size difference. Such scale differences do not change 
the qualitative rheology of dense liquids where structural 
relaxation governs the dynamics |28) . The shear rate 7 is 
switched off at i = after all systems have reached a well- 
defined NESS, imposing a constant strain for all t > Q. 
Times are reported relative to the scale of single-particle 
motion; related to free diffusion Dp, tq = R^/Dq « 0.3s 
(HS), 4.0ms and 3.4ms (PP), or to ballistic motion and 
the potential energy scale e, tq = \J AmR? ji. (MD). (The 
different forms of short-time motion may result in a shift 
of the relevant t/ro when comparing features of the long- 
time dynamics [351 ED]-) Curves start from the corre- 
sponding steady-state value a(t = Q) = crss(7) (separately 




FIG. 2. Steady-state flow curve (yssici) (open symbols) and 
residual stress 0-00(7) (filled) determined from experiment as 
CToo ~ cr(foo) in Fig. [1] (diamonds: HS, circles: PP) and MD 
simulation (squares). A solid (dotted) line shows the residual 
stress (flow curve) for the MCT model of Fig. [l] 

measured) , which is a nonlinear function of the shear rate 
(called the flow curve). 

Stresses in dense liquids are dominated by the struc- 
tural contribution, hence slow glassy dynamics governs 
their decay at long times. In the fluid, uit) relaxes to 
zero on the structural-relaxation time scale r. Approach- 
ing the glass transition, r grows beyond the experimen- 
tal window, and intermediate plateaus develop in the 
stress relaxation. In the MCT idealization, r diverges 
as permanent local caging of particles prevents the full 
relaxation of density fluctuations. The previously shear- 
molten ideal glass is then characterized by a non-relaxing 
persistent residual stress croo(7) = lim(_>.oo cr(i) > 0. Its 
shear-rate dependence highlights the nonequilibrium na- 
ture of the glassy state attained after imposing zero-flow 
conditions: different glasses exist with the same ther- 
modynamic control parameters, but different history- 
dependent frozen-in properties. 

In simulation and experiment, the residual-stress 
plateaus slowly decay, possibly as a result of creep [3lT - 
155] . which is precluded in the present MCT calculations. 
Imposing zero-stress instead of zero-strain-rate condi- 
tions on the other hand will allow for rejuvenation effects 
that may eliminate residual stresses [32 [35] . 

The residual stress Coo (7) as a function of the pre-shear 
rate is shown in Fig. [2] (filled symbols). In experiment and 
simulation, comparable values have been determined at 
a suitable intermediate time too (marked in Fig. [T]). Also 
shown (open symbols) are some of the flow curves a^siff)- 
They exhibit two features typical for glass-forming flu- 
ids: a dynamical yield stress, Oy = lim.y^o o'ss(7) > 0, 
and a monotonic increase with increasing 7. The resid- 
ual stresses 0^0 (7) , on the other hand, typically decrease 
with increasing 7 (leading to a crossing of cr (t)-versus-t 
curves at short t/ro, see Fig. [l^): the stronger the past 
fiuidization (and hence, structural distortions [18]) of the 
glass, the more effective the stress relaxation. 

For 7 — 0, the residual stress (Too approaches ay. In 



3 




10"^ 10"' lo' lo' lo' 10"" 10"' lo' 10^ 10^ yf 



FIG. 3. As Fig. [T] but a{t)/ass as a function of jt. Glass 
(liquid) states are shown in red (blue). 



ITT-MCT, the two quantities coincide in this limit, im- 
plying that for arbitrarily slow flow, the glass attains a 
certain stress that can never relax, even after the pertur- 
bation is removed. To understand this, recall that the 
shear stress is given by a nonlinear Green-Kubo relation, 
(^xyit) — J_^'^{t')G{t,t' jl-j]) dt\ where the generalized 
shear modulus G{t,t' , [j]) is a transient correlation func- 
tion that is formed with the equilibrium ensemble aver- 
age and is only affected by external perturbations active 
between its two time arguments t' < t In the ab- 

sence of other relaxation mechanisms, the same physical 
process, shear-induced breaking of cages on a time scale 
Tj ~ 1/7, dominates the history integrals determining ay 
and CToo, and Coo ~ (Ty > results. The ISHSS model of 
ITT-MCT predicts a slight increase of Coo with increas- 
ing 7, which is only seen in the HS experiment well in 
the glass. 

The stress relaxations reveal remarkable scaling behav- 
ior close to the glass transition. Normalizing stresses by 
their flow-curve value, and rescaling time with the ini- 
tial shear rate, the relaxation curves fall into two classes 
as shown in Fig. |3] Two distinct decay patterns - one 
for the liquid, one for the glass - emerge that provide 
a clear indicator to locate the glass transition through 
a series of shear-cessation experiments. The dependence 
on "pseudo-strain" jt, as if flow persisted, is remarkable 
since the equations of motion at i > contain no refer- 
ence to the past perturbation. The scaling hence can be 
thought of as the slow decay of the NESS-contributions 
to the distribution function over which dynamical quan- 
tities are averaged. In contrast, the final decay of a{t) to 
zero in the liquid does not scale with 7, as it is governed 
by the equilibrium relaxation time. 

The SGR predicts the observed scaling with jt based 
on aging phenomena [S]. Specifically, one obtains an 
asymptotic power-law decay to zero, a{t)/ass ~ (7^)"^ 
for 7i ^ 1, where a; is a temperature- like parameter. A 
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FIG. 4. Mean-squared displacement Sr^{t,tm) starting at 
a waiting time after cessation of steady shear, in units of 
particle diameters, (a) MD simulation (large particles): for 
T = 0.14, 7t„ = and 0.675 for -yro = 3 x 10"^ together 
with the NESS and equilibrium result (violet). For T = 0.1, 
tw — curves for jro = 3 x 10~^ and 3 x 10~^ (red; right 
to left), (b) Confocal microscopy on hard-sphere PMMA col- 
loids, a,t ip — 0.56. Curves for NESS and equilibrium together 
with t„ = for -yro = 0.0038/s and 0.0076/s (Peo ~ 0.0097 
and 0.019). For the smaller shear rate, = 4s is also shown. 
The thick red line is a theory fit to this curve, (c) Schematic 
ITT-MCT model, see text, for two state points e = ±10"* 
(glass/liquid) with 7ro = 10~^, and t^, = and 'yt^ — 0.01. 
(d) Intermediate-plateau values after switch-off, , as a func- 
tion of 7tu,. A dashed line marks the equilibrium plateau 
height from MD. 

finite residual stress cr 00 is not predicted by the SGR. 

To elucidate the microscopic mechanisms at play dur- 
ing stress relaxation, we turn to the mean-squared dis- 
placements (MSD) of individual particles. This is a two- 
time average, 6r'^{t,tw) = {{r{t) — r{tw))'^)t^, governed 
by the equilibrium dynamics but averaged with respect 
to the statistical ensemble at a waiting time t^, > after 
cessation of shear. MSD in the vorticity direction, for 
various shear rates and state points, and for various i^,, 
from confocal-microscopy experiments and MD simula- 
tions are shown in Fig.|4] The t^-independent MSD mea- 
sured in the equilibrium fluid and in the sheared NESS 
are shown as a reference. They exhibit an intermediate- 
time plateau indicating transient caging of particles in 
their nearest-neighbor shells. Long-time diffusion is en- 
hanced in the NESS, a non- linear response effect that is 
the microscopic analogue of shear thinning [31] . Measur- 
ing the MSD with reference to the configuration at the 
cessation point, ~ 0, the dynamics still follows that of 
the NESS, up to a time given by jt « 0.1. This connects 
to a typical strain that causes cage-breaking [19]. The 
MSD then slowly crosses over to the equilibrium dynam- 
ics, indicating an intermediate subdiffusive plateau that 
is larger than the one connected to quiescent caging. At 
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times t — t^ ^ T, all curves become independent of i^, and 
collapse on the equilibrium curves. For MD simulations 
in the glassy state the intermediate plateau extends past 
the simulated time window. The height of this plateau 
decreases with increasing 7. 

ITT-MCT (extended to describe waiting-time de- 
pendent two-time averages for small t^, [TH 
qualitatively rationalizes these MSD. Figure |4] includes 
schematic-model results for both a glassy and a liquid 
state. A well-defined second plateau emerges in the cross- 
over from the NESS to equilibrium curves and becomes 
permanent in the glass. It can be interpreted as a sec- 
ond length scale beyond the quiescent localization length, 
arising from the competition of shear-induced fluidization 
and arrest after cessation. Intuitively, it is a "Brcms- 
weg" (stopping distance) for individual particles caused 
by the progression of the distribution function from the 
perturbed state at to the quiescent state. But note 
that inertial effects play no role in the dynamics and are, 
by construction, absent from the theory. 

The ITT-MCT model approximates the i^-dependent 
MSD as a combination of the equilibrium and NESS ones. 
Given this information from experiments, the evolution 
of 5r^(t, t^) for small t^, can be predicted. This is demon- 
strated by the fit to the HS data shown in Fig. |4]d. 

The scenario of nonequilibrium relaxation discussed 
here characteristically differs from aging dynamics in spin 
glasses [SS] [35] or that following density quenches in hard 
spheres [57H5^ . There the structural relaxation time 
T{tw) grows with sample age, so that correlation func- 
tions and related two-time averages depend on t^^ in their 
long-time part, while the short-time relaxation for in- 
creasing sample age reveals more and more of an intrin- 
sic, fu,-independent relaxation. For the MSD, this implies 
a continuous shift of the long-time diffusive asymptote 
(where 6r^ ^ t) to longer times. In contrast, in Fig. |4] all 
intermediate-iu, curves show subdiffusive transients and 
an approach to the same tu,-independent long-time diffu- 
sion found in the quiescent system. 

In conclusion, we have studied stress decays and the 
microscopic dynamics of glass-forming liquids and shear- 
melted glasses after the cessation of steady shear flow. In 
the liquid, stresses relax to zero following long transients, 
on the time scale of the quiescent equilibrium system. (In 
rheological terms, the systems exhibit "thixotropy" .) 

Finite residual stresses remain in the glass. Their 
value and the initial evolution from the steady state to 
the nonequilibrium quiescent solid, are governed by the 
pre-shear rate through long-lived memory effects. These 
memory effects cause the appearance of a microscopic 
supra-caging length scale in the waiting-time dependent 
mean-squared displacements. 

Residual stresses imply that glass is not simply char- 
acterized by its thermodynamic control variables. Differ- 
ent preparation histories result in glasses that differ sub- 
tly in their structure, and possibly also in their response 



behavior, for example their elastic moduli |3]. Within 
ITT-MCT, a pre-strain dependence of the shear modu- 
lus has been studied in the flowing steady state [40]. It 
would be enlightening to compare this to a history de- 
pendence that arises in the case of temperature quenches, 
as recently studied in computer simulations [H]. 

The theory describes history-dependent glass states 
through retarded-friction contributions to the dynamics 
of density fluctuations that are modified by the past flow. 
This goes beyond traditional near-equilibrium glass tran- 
sition theories (such as standard MCT), where the relax- 
ation of small initial perturbations induced by an exter- 
nal field is related to a Kubo correlator of equilibrium 
fluctuations. This connection is based on Onsager's re- 
gression hypothesis, and it holds if the initial perturba- 
tion obeys equilibrium linear-response theory. It is vio- 
lated in the glass (unless activated processes restore er- 
godicity) [32] , since infinitesimally small rates shear-melt 
the glass [T3]. These states relax only due to external 
driving, which is the common cause of both a dynamical 
yield stress Gy = lim.y_j.o o'ss(7) > (contradicting lin- 
ear response), and of a residual stress CToo (contradicting 
equilibrium) . 

The qualitative agreement among different systems 
and methods suggests that our discussion applies gener- 
ally to glass formers where excluded volume effects dom- 
inate, causing stresses to be given by the entropic scale 
kT/R^ . This includes dense metallic glass formers, where 
an understanding of frozen-in residual stresses is impor- 
tant for understanding material stability. 

In such systems, possibly also for colloidal systems 
deep in the glass, additional effects may need to 
be accounted for, connected to spatial and temporal 
flow heterogeneities. Then, concepts such as shear- 
transformation zones (STZ) 03] or stress avalanches 
[311I3S] (or, in the case of granular materials, force chains 
[46) ) can become important. For the data presented 
above, we have checked that no such inhomogeneities are 
detectable either in confocal microscopy or during the 
simulation runs. 
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